COMPLEX HADAMARD MATRICES OF ORDER 6: 
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CNI , Abstract. In this paper we construct a new, previously unknown four-parameter family 

(3J^' of complex Hadamard matrices of order 6, the entries of which are described by algebraic 

functions of roots of various sextic polynomials. We conjecture that the new, generic family 
' Gg 4 "* together with Karlsson's degenerate family and Tao's spectral matrix form an 

exhaustive list of complex Hadamard matrices of order 6. Such a complete characterization 
might finally lead to the solution of the famous MUB-6 problem. 
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1. Introduction 

Complex Hadamard matrices form an important family of orthogonal arrays with the 
additional unimodularity constraint imposed on their entries. These matrices obey the alge- 
braic identity HH* = nl n where * stands for the Hermitean transpose, and is the identity 
matrix of order n. They appear in various branches of mathematics frequently, including 
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linear algebra [B], coding- and operator theory [TJ [T3] and harmonic analysis [T3"l \W\ . They 
play an important role in quantum optics, high-energy physics, and they are one of the key 
ingredients to quantum teleportation- and dense coding schemes [20] and mutually unbiased 
bases (MUBs) |2J. An example of complex Hadamard matrices is the Fourier matrix F n , 
well-known to exists for all n. It is natural to ask how does a "typical" complex Hadamard 
matrix of order n look like, and a satisfying answer to this question can be given provided 
we have a complete characterization of Hadamard matrices of order n at our disposal. These 
types of problems, however, are notoriously difficult even for small n. Naturally, one is inter- 
ested in the essentially different matrices only, and we identify two matrices H and K and 
say that they are equivalent if H = P1D1KD2P2 for some unitary diagonal matrices Di, D 2 
and permutational matrices P\,P%. Recall that a complex Hadamard matrix is dephased, 
if all entries in its first row and column are equal to 1. While studying and classifying real 
Hadamard matrices is naturally a discrete, finite problem which can be handled by deep 
algebraic methods and sophisticated computer programs to some extent, the complex case, 
however, behaves essentially different. In particular, due to the appearance of various para- 
metric families one cannot hope for a finite list of inequivalent matrices, but rather for a 
finite list of constructions, each of them leading to an infinite family of complex Hadamard 
matrices. 

The complete classification of complex Hadamard matrices is available up to order n = 5 
only. It is trivial that F n is the unique complex Hadamard matrix for orders n < 3. The case 



Date: August, 2010., Preprint. 

This work was supported by the Hungarian National Research Fund OTKA K-77748. 

1 



2 FERENC SZOLLOSI 

n = 4 is still elementary, and it was shown by Craigen that all complex Hadamard matrices 
of order 4 belong to an infinite, continuous one-parameter family [4]. In order 5 we have 
uniqueness again, a result which is absolutely non-trivial already. In particular, Lovasz was 
the first who showed [12] that F§ is the only circulant complex Hadamard matrix in this 
order, and a decade later Haagerup managed to prove the uniqueness of F 5 by discovering 
an algebraic identity (cf. formula ( 1171) ) relating the matrix entries in a surprising way [7]. 

In order 6 various one- [TJ [HI [21] , two- [TH [T7] and three-parameter families [91 [10] have 
been constructed recently and it is conjectured that these are part of a more general, four- 
parameter family of complex Hadamard matrices, yet to be discovered [2]. This conjecture 
is supported by overwhelming numerical evidence [16], however so far only a fairly small 
subset of it was described by closed analytic formulae, including an isolated matrix Sq and 
a three-parameter matrix K\ [9j [10] . 

The reason why the 6x6 case received significant attention in the past couple of years 
is the fact that complex Hadamard matrices are closely related to mutually unbiased bases. 
Recall that two orthonormal bases of C d , B\ and B2 are unbiased if for every e £ B%, f G B2 
we have |(e, f)\ 2 = 1/d. A family of orthonormal bases is said to be mutually unbiased if 
every two of them are unbiased. The famous MUB-6 problem asks for the maximal number 
of mutually unbiased bases in C 6 . On the one hand this number is at least 3, as there 
exists various infinite families of triplets of MUBs in this order [HI [TTJ [21], on the other 
hand it is well-known that it cannot be larger than 7 (cf. the references of [2]). In fact, 
it is conjectured that a triplet is the best one can come up with in dimension 6 |21j . The 
connection between MUBs and Hadamard matrices of order 6 has been exploited in [S] very 
recently, where a discretization scheme was offered to attack the problem and it was proved 

by means of computers, however, in a mathematically rigorous way, that the members of the 

(2) 

two-parameter Fourier family F 6 (a, b) and its transpose cannot belong to a configuration of 
7 MUBs containing the standard basis in dimension 6. One reasonable hope to finally settle 
the MUB-6 problem is to give a complete characterization of complex Hadamard matrices 
of order 6 and apply the same technique to them. 

The goal of this paper is to propose a general framework towards the complete classification 
of complex Hadamard matrices of order 6. In particular, by characterizing the orthogonal 
triplets of rows in complex Hadamard matrices we generalize an observation of Haagerup [7] 
to obtain a new algebraic identity relating the matrix entries in an unexpected way. This is 
an essentially new tool to study complex Hadamard matrices of small orders, and one of the 
main achievements of this paper. We apply this result to obtain complex Hadamard matrices, 
moreover we conjecture that the the construction we present here reflects the true nature 
of complex Hadamard matrices of order 6. It has the following three features: Firstly, it is 
general in contrast with the earlier attempts where always some additional extra structure 
was imposed on the matrices including self-adjointness p], symmetry [14] , circulant block 
structure [TF] or if 2 -reducibility [;9j . Secondly, it has 4 degrees of freedom and thirdly all the 
entries of the obtained matrices can be described by algebraic functions of roots of various 
sextic polynomials. This suggests on the one hand the existence of a four-parameter family 
of complex Hadamard matrices of order 6 and reminds us on the other hand the fact that 
the desired algebraical description where the entries are expressed by radicals might not be 
possible at all. However, from the applicational point of view, and in particular, to utilize the 
computer-aided attack of [8] to the MUB-6 problem we shall need these matrices numerically 
anyway. 
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The outline of the paper is as follows. In Section 2 we briefly discuss the main ideas 
of the construction to motivate the various auxiliary results we prove there. The excited 
reader might want to skip this section at first, and jump right ahead to Section 3, where the 
construction of the new family is presented from a high-level perspective. In Section 4 we 
analyze the construction thoroughly. 



2. Preliminary results 



In this section we present the ingredients necessary to construct our new family of complex 
Hadamard matrices of order 6, including a characterization of the mutual orthogonality of 
three rows in Hadamard matrices (cf. Theorem l2.4p . First, however, we would like to motivate 
these efforts by describing the main ideas of the construction. 

We start with a submatrix 



(1) 



E(a, b, c, d) :- 



1 1 1 
lab 
led 



and attempt to embed it into a complex Hadamard matrix of order 6 
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with 3x3 blocks E, B, C and D in two steps, as follows. First we construct the submatrices 
B and C featuring unimodular entries to obtain three orthogonal rows and columns of G$. 
Secondly we find the unique lower right submatrix D to get a unitary matrix. Should the 
entries of this matrix become unimodular, we have found a complex Hadamard matrix. We 
conjecture that the submatrix E can be chosen, up to equivalence, in a way that there will 
be only finitely many candidates for the blocks B and C and therefore we can ultimately 
decide whether the submatrix E can be embedded into a complex Hadamard matrix. The 
resulting matrix G§ can be thought as the "Hadamard dilation" of the operator E. 

We shall heavily use the following through the paper without any further comment: the 
conjugate of a complex number of modulus 1 is its reciprocal, and hence the conjugate of 
a multivariate polynomial with real coefficients depending on indeterminates of modulus 1 
is just the polynomial formed by entrywise reciprocal of the aforementioned indeterminates. 
We computed various Grobner bases [3] in this paper with the aid of Mathematica. The 
reader is advised to use a computer algebra system for bookkeeping purposes and consult 

(0) ^(3) 



for the standard notations for well-known complex Hadamard matrices such as S$ , K { 



etc. 



We begin with recalling two elementary results from the existing literature. 



Lemma 2.1. Suppose that we have a partial row (1, a, b, e, *, *) composed from unimodular 
entries. Then one can specify some unimodular numbers s\ and s 2 in place of the unknown 
numbers * to make this row orthogonal to (1, 1, 1, 1, 1, 1) if and only if 



(3) 



II + a + 6 + el < 2. 



4 FERENC SZOLLOSI 

Proof. To ensure orthogonality, we need to have l + a + b + e + si + S2 = from which it 
follows that |1 + a + b + e\ = |si + s 2 | < 2. It is easily seen geometrically, that in this case 
we can define the unimodular numbers required. □ 



The missing coordinates featuring in Lemma [2.14 Si and S2, can be obtained algebraically 
through the well-known 

Lemma 2.2 ((Decomposition formula, [14])). Suppose that the rows (1,1,1,1,1,1) and 
(1, a, b, e, si, S2) containing unimodular entries are orthogonal. Let us denote by S := 1 + a + 
6 + e, and suppose that < |S| < 2. Then 



(4) ^=-| ±i iljV 1 - 

//£ = £/ien Si zs independent from a,b,e but S2 = — si- 

Proof. Clearly Si and S2 are the unimodular numbers with s± + S2 = — S. □ 

Now we proceed by investigating the orthogonality of triplets of rows. In order to do this, 
the following is a crucial 

Definition 2.3 ((Haagerup polynomial)). The Haagerup polynomial % associated to the 
rows (1, 1, 1, 1, 1, 1), (1, a, b, e, *, *) and (1, c, d, f, *, *) of a complex Hadamard matrix read 

H(a,b,c,d,e,f) := (1 + a + b + e)(l + c + d + J)(l + ca + db + ft). 

The first result of ours is the following 

Theorem 2.4. Suppose that we have the partial rows (1, a, b, e, *, *) and (1, c, d, /, *, *), 
composed from unimodular entries. Then one can specify some unimodular numbers si, s 2 , S3 
and S4 in place of the unknown numbers * to make these rows together with (1, 1, 1, 1, 1, 1) 
mutually orthogonal if and only if 

(5) K(a,b,c } d,e } f) = 4 - |1 + a + b + e\ 2 - |1 + c + d + f\ 2 - |1 + ca + db + fe\ 2 
with 

(6) \n(a,b,c,d,e,f)\<8. 

Proof. First we start by proving that (J5]) holds. To do this, we utilize Haagerup's idea j7] as 
follows: by pairwise orthogonality, we find that 

l + a + b + e= —si — s 2 
1 + c + d+ f = —s 3 - s 4 
1 + ca + db + fe = — s 3 s"i — s 4 s 2 
Now, by multiplying these three equations together we find that ([6]) follows and further 

H = -(Si + S 2 )(s 3 + S 4 )(s 3 Si + S4S2) = -(S1S3 + S 2 S 4 + S1S4 + S2S 3 )(s 3 Si + S4S2) = 

-|s!S 3 + s 2 s 4 | 2 - (sis 4 + s 2 s 3 )(s 3 si + S4S2) = -|sis 3 + s 2 s 4 | 2 - 2K(sis 2 + s 3 s 4 ). 
To conclude the proof, we need to show that 

2K(sis 2 + s 3 s 4 ) = |l + a + 6 + e| 2 + |l + c + rf + /| 2 -4 
holds, however this follows from the Decomposition formula easily. 
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To see the converse direction, we need to show that (jSJ) essentially encodes orthogonality. 
Let us use the notations E := 1 + a + b + e, A := 1 + c + d + f, * := 1 + ca + db + fe. With 
this notation condition (JHJ) boils down to 

(7) K = SA* = 4 - |E| 2 - |A| 2 - |*| 2 . 

Clearly, if |E| < 2 and |A| < 2 hold, then by the decomposition formula we can find si,S2j 
s 3 and s 4 to ensure orthogonality to row (1,1,1,1,1,1). Now observe, that the mutual 
orthogonality of rows (1, a, b, e, Si, s 2 ) and (1, c, d, f, s 3 , s 4 ) reads 

(8) * + s 3 S! + s 4 s 2 = 0. 

Suppose first that we have the trivial case E = A = 0. Then, by the decomposition formula 
we have s 2 = — si and s 4 = —S3, and ([7j) implies that = 2. Therefore, if we set the 
unimodular number S3 := —^51/2 the orthogonality equation (jSJ) is fulfilled. 

Suppose secondly, that we have A = 0, but S ^ 0. Then we have s 4 = — s 3 , and from (j7j) 
it follows that |E| < 2, and in particular 




(9) |*| = V4- |S| 2 . 

Now we can use the Decomposition formula to find out the values of s\ and s 2 and the 
orthogonality equation (JS} becomes 

(10) * + S3 ^_ 2i ^^/l-^j =0 

This holds, independently of S3, if |E| = 2, as by * = follows. Otherwise, set the 
unimodular number 



to ensure the orthogonality through (llOj) . 

Finally, let us suppose that E 7^ and A^0. Now observe that in this case the value of 
* needed for formula © can be calculated through (JTj). The other ingredient, namely the 
value of S3S1 + s 4 s 2 can be established through the Decomposition formula, once we derive 
the required bounds |E| < 2 and |A| < 2. Depending on the value of "H, we treat several 
cases differently. 

CASE 1: Suppose that -|*| 2 < H. This implies, by formula ©, that |E| 2 + |A| 2 < 4, 
and in particular |E| < 2, |A| < 2 hold. Next we calculate |*| from (^). 
Suppose first that % > 0. Hence, after taking absolute values, ([7j) becomes 

|*| 2 + |E||A||*| + |E| 2 + |A| 2 -4 = 0, 

and the only non-negative root is 

-|E||A| + y/(A- |E| 2 )(4- |A| 2 ) 



(11) 1*1 = 
v ) 11 2 

Now suppose that — 1*| 2 < Tl < 0. Hence, after taking absolute values, (J7|) becomes 

(12) |*| 2 - |E||A||*| + |E| 2 + |A| 2 - 4 = 0, 

and we find that the only non- negative root we have under the assumption |E| 2 + |A| 2 < 4 
reads 

IEIIAI + J(A - IEPV4- IAR 
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CASE 2: Suppose now that H < -|\I/| 2 . This implies that |E| 2 + |A| 2 > 4, and we do not 
have a priori the bounds |E| < 2, |A| < 2. Nevertheless, we derive equation ffT2l again, and 
we find that the values of \$!\ can be any of 



(14) \m 



1,2 



|E||A| ± ^(4- |E| 2 )(4- |A| 



provided that the roots are real, namely we have either |E| > 2 and |A| > 2 or |E| < 2 and 
|A| < 2. The first case is, however, not possible, as it would imply \^\ > 2 contradicting the 
crucial condition ([6]). 

Once we have established the bounds |S| < 2, |A| < 2 and the value(s) of \^\ has been 
found, we are free to use the Decomposition formula to obtain the values of s\, s 2 , S3 and 
s 4 . Clearly, we can set Si with the + sign, while s 2 with the — sign as in formula (0J, up to 
equivalence. However, we do not know a priori how to distribute the signs amongst S3 and 
S4, and to simplify the notations we define 



A . A / |A| 2 A . A / |A| 2 

(15) S3 = ±1— -Wl — J -, S 4 = T I7-T-7A/ 1 — - — -. 

v ; 3 2 |A| V 4 ' 2 |A| V 4 

In particular, by using ([7]) and (Tl5]) we find that the orthogonality equation ([8]) becomes 



, , 4 — E 2 — A 2 - * 2 EA EA / E 2 / A 2 

16 — — + ±2— — T 4/i-J-J-Wi-L_L = o, 

y J EA 2 |E||A|V 4 V 4 

where the ± sign agrees with the definition of S3. To conclude the theorem plug in all of the 
possible values of as described in fTTTT) . ffT3l) . ffT4j) into ffl~6l) to verify that for some choice 
of the sign it holds identically. □ 

Remark 2.5. The two possible signs described by formula ( Fl4"l) can be realized. In particular, 
there are two different orthogonal triplet of rows composed of sixth roots of unity where 
|E| = |A| = in both cases, however, in one of the cases = 1 while = 2 in the 
other. 

Corollary 2.6 ((Haagerup's trick, [7J)). Suppose that the rows (1, 1, 1, 1, 1, 1), (1, a, b, e, Si, s 2 ) 

and (1, c, d, f, s 3 , s 4 ) composed of unimodular entries are mutually orthogonal. Then 

(17) H(a,b,c,d,e,f) G K. 



Haagerup used the property ( 1171) to give a complete characterization of complex Hadamard 
matrices of order 5, or equivalently, describe the orthogonal maximal abelian *-subalgebras 
of the 5x5 matrices [7]. Since then it was used in [T] and [TJ] to construct new, previ- 
ously unknown complex Hadamard matrices of order 6 as well. However, to guarantee the 
mutual orthogonality of three rows the necessary condition ( ITTj) should be replaced by the 
more stronger identity ([5]). Nevertheless, ( TT7j) will play an essential role in this paper too. 
These type of identities are extremely useful as they feature less variables than the standard 
orthogonality equations considerably simplifying the calculations required. 

Solving the system of equations (15|)-( TT7|) is the key step to obtain the submatrices B and 
C of Gq. Once we have three orthogonal rows and columns we readily fill out the remaining 
lower right submatrix D. This is explained by the following two lemmata, the first of which 
being a special case of a more general matrix inversion 
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Lemma 2.7. If U and V are n x n matrices then 

(i n + uvy 1 = i n -u (/ n + vuy 1 v 

provided that one of the matrices I n + UV or I n + VU is nonsingular. 

Proof. By symmetry, we can suppose that the matrix I n + VU is nonsingular. Then, we have 

(/„ + uv)(i n - u(i n + vuy l v) = i n + uv- u(i n + vu)(i n + vuy l v = i n . 

□ 

Lemma 2.8. Suppose that we have a 6 x 6 partial complex Hadamard matrix consisting of 
three orthogonal rows and columns, containing no vanishing 3x3 minor. Then there is a 
unique way to construct a unitary matrix containing these rows and columns as a submatrix. 

Proof. Let U be a 6 x 6 matrix with 3x3 blocks A, B, C and D, as the following: 

A B 



U 



C D 



By the orthogonality of the first three rows and columns and using the fact that the entries 
are unimodular, we have 

(18) AA* + BB* = 6/ 3 

(19) A* A + C*C = 6/3 

To ensure orthogonality in-between the first three and the last three rows we need to have 
AC* + BD* = O3, the all matrix. As B is nonsingular by our assumptions we can define 

(20) D := -CA*(B- 1 )*. 

Now we need to show that the last three rows are mutually orthogonal as well. Indeed, by 
using (!20l) and (jTBl we have 

CC* + DD* = C{h + A*(BB*y 1 A) C* = C (/ 3 + A* (6/ 3 - AA*)' 1 A) C*, 

which, by Lemma [2.71 and (TT9T) is 




1 



-1 



= c ( h - q a * a ) c * = QC (c*cy L c* = 6/3. 

□ 

We do not state that the obtained unitary matrix U is Hadamard, which is not true in 
general. Recall that our goal is to embed the submatrix E into the matrix G 6 (cf. ([I])-®). 
We have the following trivial 

Lemma 2.9. Suppose that a submatrix E can be embedded into a complex Hadamard matrix 
Gq of order 6 in which the upper right submatrix B is invertible. Then for some unimodular 
submatrix C for which the first three columns of G 6 are orthogonal the lower right submatrix 
D = —CE*{B~ 1 )* is unimodular. 

Proof. Indeed, this is exactly what embedding means. □ 

In particular, if the submatrices B and C are chosen carefully, the unimodular property 
of D follows for free. 
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Corollary 2.10. Start from a submatrix E and suppose that there are only finitely many 
(invertible) candidate matrices B G SOL B and C G SOL c such that the first three rows and 
columns of the matrix G e are orthogonal. Then E can be embedded into a complex Hadamard 
matrix of order 6 if and only if there is some B G SO Lb and C G SOLc such that the matrix 
D = —CE*(B~ 1 )* is unimodular. 

Note that due to the finiteness condition in Corollary 12.101 once we have all (but finitely 
many) candidate matrices B and C we can decide algorithmically whether the submatrix E 
can be embedded into a complex Hadamard matrix. 

The next step is to characterize 6x6 complex Hadamard matrices with vanishing 3x3 
minors. To do this we need two auxiliary results first. 

Lemma 2.11. Suppose that in a dephased 6x6 complex Hadamard matrix there exist a 
noninitial row (or column) containing three identical entries x. Then x = ±1 and this row 
(or column) reads (1, 1, 1, —1, —1, —1), up to permutations. 

Proof. Suppose to the contrary, that there are three nonreal numbers x in a row (column). 
Then the sum of these numbers x together with the leading 1 read |l + 3x| 2 = 10 + 63ft(x) > 4, 
and hence, by Lemma \2. II this row (column) cannot be orthogonal to the first row (column), 
a contradiction. To ensure orthogonality to the first row (column) we should specify the 
remaining three entries to — x and hence the last part of the statement follows. □ 

Recall that the core of a dephased complex Hadamard matrix of order n is its lower right 
(n — 1) x (n — 1) submatrix. A vanishing sum of order k is a fc-term sum adding up to 0. 
The following breakthrough result was obtained very recently. 

Theorem 2.12 ((Karlsson, [9J HO]))- Let H be a dephased complex Hadamard matrix of 
order 6. Then the following are equivalent: 

(a) H belongs to the three-parameter degenerate family Kq 3 \ 

(b) H is H 2 -reducible; 

(c) The core of H contains a — 1; 

(d) Some row or column of H contains a vanishing sum of order 2; 

(e) Some row or column of H contains a vanishing sum of order 4. 

In particular, Theorem 12.121 gives a characterization of complex Hadamard matrices of 
order 6 containing F 2 as a submatrix. The term if 2 - re ducibility refers to the beautiful 
structure of these matrices: they have a canonical form in which all 9 of their 2x2 submatrices 
are complex Hadamard. Part (c) and (d) of Theorem 12.121 allow us to quickly recognize if 

(3) 

a matrix belongs to the family K$ , and we shall heavily use these conditions through our 
paper. Part (e) tells us that once we have four entries in a row or column of a matrix which 

(3) 

lies outside the family K & the Decomposition formula readily derives the unique remaining 
two values through (J4]). As the family forms a three-parameter subset, the matrices it 
contains are atypical, hence the adjective "degenerate" . 

Corollary 2.13. Suppose that in a dephased 6x6 complex Hadamard matrix H there exist a 
noninitial row (or column) containing three identical entries. Then H belongs to the family 



Lemma 2.14. Suppose that a 6 x 6 complex Hadamard matrix H has a vanishing 3x3 
minor. Then H belongs to the family . 
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Proof. Suppose that H has a vanishing 3x3 minor, say the upper left submatrix E(a, b, c, d), 
as in formula (pQ). Such an assumption can be made, up to equivalence. As det(E) = 
b + c — a — d + ad — be = 0, we find that if any of the indeterminates a,b,c,d is equal to 1 
then E contains a noninitial row (or column) containing three Is, and therefore by Corollary 
12.131 we conclude that this matrix belongs to the family K^p . Otherwise, we can suppose 
that none of a, b, c, d is equal to 1 and hence we find that d — (a + be — b — c)/(a — 1), which 
should be of modulus one. To ensure this solve the equation dd — 1 = to find that either 
b = a or c = a should hold, but then we have either d = c or d = b as well. In particular, we 
find that E has two identical rows (or columns). After enphasing the matrix, again, we find 
that there is a full column (or row) of entries 1 and a reference to Corollary 12.131 concludes 
the lemma. □ 

(3) 

Therefore to investigate those matrices which lie outside the family K§ we can safely use 
Lemma 12.81 and in particular the inversion formula (l2~Uj) . 

It turns out, that the isolated matrix (cf. [IB]) requires a special treatment as well. It 
is featured in the following 

Lemma 2.15. Suppose that in a 6 x 6 dephased complex Hadaraard matrix H there is a 
noninitial row and column composed of cubic roots of unity. Then H is either equivalent to 



Sf^ or belongs to the family K§ 

Let us denote by uj the principal cubic root of unity once and for all, that is uj := e 



(3) 



2yri/3 



Proof. First suppose that the cubic row and column meet in a common 1. Then our matrix 
looks like as the matrix H on the left below, up to equivalence: 





■ 1 


1 1 1 


1 


1 ■ 




■ 1 


1 1 


1 


1 


1 ■ 
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1 UJ UJ 


UJ 2 






1 


1 u 


UJ 


UJ 2 


UJ 2 


H = 


1 
1 


uj a b 
uj * * 


c 
* 


d 

* 


, H' = 


1 
1 


a 1 

* u 


6 

* 


c 
* 


d 

* 




1 


UJ 2 * * 


* 


* 




1 


* UJ 2 


* 


* 


* 




. 1 


uj 2 * * 


* 


* 




. 1 


* UJ 2 


* 


* 


* 


Now by orthogonality of the first three 


rows we find that 


a + b = 


— UJ 


and 


c + d 


we can assume, up to equivalence 


, that a 


= 1,6 = - 




and c = 


= UJ 


, d - 


-- UJ 2 . 



But then, we 

can fill out the fourth row, and the third and fourth column as well. We conclude that the 
obtained matrix is equivalent to Sf\ 

Secondly, let us suppose that the cubic row and column meet in a common uj. This matrix 
H' is depicted on the right above. Then, by calculating the orthogonality equations, we find 
that either a = — 1, and hence the matrix (if it can be completed to a Hadamard at all) 



This 



belongs to the family K 6 , or a = uj, b = uj 2 , and, up to equivalence c = uj, d = uj y 
implies that the third row and third column feature cubic entries only meeting in a common 
1 therefore reducing the situation to the first case. The third case, namely when the cubic 
row and column meet in a common uj 2 can be treated similarly. □ 

Now we turn to the presetting of the submatrix E(a,b,c,d) (see flTJ). In order to avoid 
the case when the system of equations (151)- (1171) is linearly dependent we need to exclude 
various input quadruples (a,b,c,d). However, it shall turn out that we are free to do such 
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restrictions, up to equivalence. Before proceeding further, let us define the following two- 
variable function mapping T 2 to C as follows: 

£(x, y) := x + y + x 2 + y 2 + xy 2 + x 2 y. 

We say that y is an elliptical pair of x, if 8(x,y) = 0. Observe that for a given i / -1 
the sum of its elliptical pairs read yi + y-z = — (1 + x 2 )/(l + x). The following is a strictly 
technical 

Proposition 2.16 ((Canonical transformation)). Suppose that we have a complex Hadamard 
matrix H inequivalent from and any of the members of the family . Then H has a 
3x3 submatrix E(a,b,c,d) as in formula ([!]), up to equivalence, satisfying 

(21) (6 - l)(c - 1)(6 - d 2 ){c - d 2 )(b - c){bc - d)S{b, d)£{c, d) ^ 0. 

Proof. The strategy of the proof is the following: first we pick a "central element" d from 
the core of the matrix and then we show that b and c can be set satisfying (j2"T|) . Recall, that 
by Lemma 12.151 there is no a noninitial row and column composed from cubic roots of unity 
in the matrix. 

First let us assume that there is a 1 in the core somewhere. Suppose that there is a, to 
say, row full of cubics. In this case set d = uj, and c = uj 2 . Now in the column containing d 
there is a non-cubic entry 7, and we are free to set 6 = 7. If there is neither a full row nor a 
column of cubics in the matrix, then set d — 1, and choose a non-cubic c from its row. Note 
that there cannot be a further noninitial 1 in the row or column of d by Corollary 12.131 Now 
observe that as the elliptical pairs of 1 are uj and uj 2 , we can choose a suitable b from the 
column of d unless all entries there are members of the set {uj, uj 2 , c, c}. Note that uj together 
with uj 2 cannot be in the column of d at the same time, and from this it is easily seen that 
we cannot define the value of b only if the column of d is one of the following four cases, up 
to permutations: (1, 1, cu, c, c, c), (1, 1, uj, c, c, c), (1, 1, uj 2 , c, c, c), (1, 1, uj 2 , c, c, c). However, 
by normalization and by orthogonality, the sum of the entries in this column should add 
up to 0, and we find in all cases that the unimodular solution to c is a cubic root of unity, 
contradicting the choice of c. Therefore one of the entries in the column of d is different from 
uj, uj 2 , c,c which will be chosen as b. 

Secondly, let us suppose, that there is no 1 in the core. In particular, all entries in the 
core are different rowwise and columnwise. 

Pick any d from the core of the matrix. Let us denote by c\ and C2 the elliptical pairs 
of d (maybe c\ = C2, or they are undefined). Now we have several cases depending on the 
appearance of these values in the row and column containing d. 

CASE 1: C\ and C2 present in both the row and column containing d. Hence, in the row 
and column containing d there are the entries 1, d, C\ and C2 already, the remaining two, 
a and f3 are uniquely determined by the Decomposition formula. Note that a 7^ d 2 , as 
otherwise we would have /3 = —(2d + d 2 + d 3 ) /(1 + d) (as the sum of all entries in a noninitial 
row add up to 0, and the sum of the elliptical pairs is known), which is unimodular if and 
only if d = ±i or d = uj or d = uj 2 . In the first case we have (3 = =Fi and we are dealing with 
a member of the family K & . The second and third case imply that we have a full row and a 
full column of cubics, a contradiction. If ^ d/at then by picking c = a from the row we can 
set b = f3 in the column. Otherwise, in case we have (3 = d/a, then reset the central element 
d to the a which is in the same row as d. Now observe that after this exchange we should 
met the requirements of Case 1 again (otherwise we are done here), and hence the elliptical 
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pairs of a, a\ and a 2 , should present in the row and column of a, which therefore contain 
exactly the same entries. However a 12 7^ d, as otherwise orthogonality with the condition 
£(a, d) = would imply d = ±1, a contradiction; a 12 7^ d/a, as otherwise d = u or d = u 2 
would follow from the same argument implying that the row containing a has a noninitial 
1, a contradiction. Therefore, the only option left is that a^ = c\^- But then we can set 
c = d 7^ a 2 and b = d/a 7^ a 2 , and we are done. 

CASE 2: ci and c 2 present in (to say) the row containing d, but only one of these values 
(say Ci) is present in the column of d. Let us denote by a and /3 the two further entries in 
this row which, again, are different from d 2 . In the column of d there is already 1, d and 
Ci, and observe that the remaining three entries cannot be (d 2 ,a,f3) as this would imply 
d 2 = C2, contradicting our case- assumption. Therefore one of the three unspecified entries 7 
is different from d 2 , a and (3. Now if 7 7^ d/a then set c = a, 6 = 7 otherwise set c = (3, 
6 = 7. We are done. 

CASE 3: Only the value C\ is in (to say) the row containing d and in the column of d as 

well. This is a tricky case, as it might happen that the undetermined triplet in both the 

d-th row and column is precisely (d 2 , a, d/a), and therefore we cannot ensure condition ([21]) . 

However, from the orthogonality equation 1 + d + d 2 + c\ + a + d/a = 0, from its conjugate, 

and the elliptical condition £{c\,d) = we can form a system of equations, the solution of 

which can be found by computing a Grobner basis. By investigating the results, we find 

(3) — 

that either a = —1 or a + d = leading us to the family K§ or c% = d. In the last case, 
however, we can calculate the values of d and a explicitly. In particular, we find that the 
values of d and a are given by some of the unimodular roots of the following polynomials 
l + 2d + 2d 3 + d 4 = and 1 + 4a - 2a 2 - 8a 3 - 8a 4 - 8a 5 - 2a 6 + 4a 7 + a 8 = 0. Now reset the 
"central" entry d to a, and observe that the elliptical pairs of a are not present in its row, 
and therefore the conditions of this subcase are no longer met. Otherwise we can suppose 
that the undetermined triplet in the column of d is not (d 2 , a, d/a). Set c = a 7^ d 2 . Pick 7 
from the column which is different from d 2 , a, d/a, set 6 = 7 and we are done. 

CASE 4: Only the value C\ is in (to say) the row containing d and in the column of d 
there is the other elliptical value c 2 7^ c\. We can suppose that two of the undetermined 
entries in the row of d satisfy a 7^ d 2 and (3 7^ d 2 and set c = a. Now if in the column the 
undefined triplet is precisely (a, d/a, d 2 ), then observe that the same triplet cannot appear 
in the row, as otherwise c\ = c 2 would follow. Therefore we can reset c to a value different 
from a, d/a, d 2 , and set b = a. Otherwise there is an entry in the column which can be set 
to b, we are done. 

CASE 5: In the column of d there is no elliptical value at all. Pick any c = a 7^ d 2 from 
the row. Now in the column there are four unspecified entries. Clearly, one of them, say 7 
will be different from d 2 , a and d/a. Set 6 = 7. We are done. □ 

Not every submatrix E can be embedded into a complex Hadamard matrix of order 6. 
To offer a necessary condition, let us recall first that an operator A is called a contraction, 
if || A || 2 < 1, where ||.|| 2 denotes both the Euclidean norm on C 6 and the induced operator 
norm on the space of 6 x 6 matrices. We have the following 

Lemma 2.17. If A is any 3x3 submatrix of a complex Hadamard matrix H of order 6 then 
A/\f§ is a contraction. 
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Proof. Clearly, we can assume that this submatrix A is the upper left of the matrix H, which 
we will write in block form, as follows: 



Now suppose, to the contrary that there is some vector s, such that ||As|| 2 > v6||s|| 2 and 
consider the block vector s' := (s, 0) T G C 6 . We have 

\\Hs'\\ 2 = ||(As,Cs) T || 2 > ||(As,0) T || 2 = \\As\\ 2 > V6\\s\\ 2 = V6||s'|| 2 = \\Hs'\\ 2 , 

where in the last step we used that the matrix H/ V6 is unitary. □ 

In particular, we have the following 

Corollary 2.18. If the submatrix E can be embedded into a complex Hadamard matrix of 
order 6, then every eigenvalue A of the matrix E*E satisfy A < 6. 

Corollary 12.181 is a useful criterion to show that a matrix E cannot be embedded into 
a complex Hadamard matrix, however, it is unclear how to utilize it for our purposes. In 
particular, we do not know how to characterize those 3x3 matrices which satisfy its con- 
ditions. Also it is natural to ask whether the presence of the large eigenvalues is the only 
obstruction forbidding the submatrix E to be embedded. The answer to this question might 
depend on the dimension, as it is easily seen that while every 2x2 matrix can be embedded 
into a complex Hadamard matrix of order 4, only a handful of very special 2x2 matrices 
can be embedded into a complex Hadamard matrix of order 5 due to the finiteness result of 
Haagerup [7]. 

Now we are ready to present a new, previously unknown family of complex Hadamard 
matrices. The next section gives an overview of the results. 



3. The construction: A high-level perspective 

Here we describe the generic family G^ from a high-level perspective. In particular, we 
outline the main steps only, and do not discuss some degenerate cases which might come up 
during the construction. The next section is dedicated to investigate the process in details. 
The main result of this paper is the following 

Construction 3.1 ((The Dilation Algorithm)). Do the following step by step to obtain 
complex Hadamard matrices of order 6. 

#1: INPUT: the quadruple (a, b, c, d), forming the upper left 3x3 submatrix E(a, b, c, d), as 
in formula 

#2: Use Haagerup's trick to the first three rows of G^ (see (jSJ) to obtain a quadratic 
equation to /: 

(22) ^i+^ 2 / + ^3/ 2 = 0, 

where the coefficients J-\, Ti and J-3 depend on the parameters a, b, c, d and the inde- 
terminate e, and derive the following linearization formula from it: 
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#3: Use Theorem 12.41 to obtain another quadratic equation to /: 

(24) g l + g 2 f + g 3 f = o, 

where, again, the coefficients Qi, Q2 and Q% depend on the parameters a, b, c, d and the 
indeterminate e. Plug the linearization formula ( 1231) into ( 1241) and rearrange to obtain 
the companion value of e f — F(e), where 

(25) F(e) := ■ 

#4: As I /| = 1 should hold, calculate the sextic polynomial G(e) coming from the equation 
F(e)F(e) — 1 = and solve it for e. 

#5: Amongst the roots of G find all unimodular triplets (e, s±, S2) satisfying e + si + S2 = 
— 1 — a — b, calculate the companion values / = F(e), s 3 = F(si) and s 4 = F(s 2 ) 
through formula (1251) and store all sextuples (e, si, s 2 , /, s 3 , s 4 ) in a solution set called 

#6: Repeat steps #2-#6 to the transposed matrix (i. e. to the first three columns), mutatis 
mutandis to obtain the solution set SOLq- 

#7: For every pair of sextuples from SO Lb and SOLq construct the submatrices B and 
C, check if the first three rows and columns are mutually orthogonal and finally use 
Lemma [2.81 to compute the lower right submatrix D through formula (1201) . 

#8: OUTPUT: all unimodular matrices found in step #7. 

Construction 13.11 gives the essence of the new family discovered, and in the next section 
we shall give it a mathematically rigorous, low-level look. 



4. The construction: The nasty details 

Here we investigate the steps of Construction 13.11 in details. 

STEP #1: Choose a quadruple (a,b,c,d) in compliance with the Canonical Transformation 
described by Proposition 12.161 as the INPUT, and form the submatrix E. Check if it meets 
the requirements of Corollary 12. 18j if yes, then proceed, otherwise conclude that it cannot be 
embedded into a complex Hadamard matrix of order 6. Experimental results show that once 
three out of the four parameters are fixed the last one can be easily set to a value such that 
the quadruple (a, b, c, d) leads to a complex Hadamard matrix. Heuristically this means that 
there is nothing "mystical" in the choice of the initial quadruple and hence the parameters 
should be independent from each other. 

STEP #2: To obtain formula (J2*3"j) we need to see that J-3 ^ 0, independently of e. Indeed, 
suppose otherwise, which means that the following system of equations (where the last two 
are the conjugate of the first two, up to some irrelevant constant factors) 

abc + a 2 bc + abd + ab 2 d + a 2 bcd + ab 2 cd = 0, 

b 2 c + ab 2 c + a 2 d + a 2 bd + acd + bed = 0, 

a + b + ac + abc + bd + abd = 0, 

a 2 b + ab 2 + be + b 2 c + ad + a 2 d = 0, 

are fulfilled. We compute a Grobner basis and find that the polynomial 



bc{\ + c 2 )(c - d)d(l + d 2 ){c 2 + d 2 )(l + d + d 2 ) 
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is a member of it. After substituting back into the original equations we find that there is 
either a vanishing sum of order 2inEora = b= l and therefore the whole family is a 

(3) 

member of Kg , or we have E = F% or E = F 3 * but these matrices have b = c which however 
is not allowed by the Canonical Transformation. 

It might happen that ^ but there is a unimodular e making it vanish, which cannot 
be anything else, but 

. , a 2 b + a 2 d + ab 2 + ad + b 2 c + be 

26 e = — i n n — T~- 

abc + aba + ac + a + bd + b 

Nevertheless, we can suppose that in case of one of the pairs (e, /), (si, S3) and (S2, S4) we do 
not set e as above, otherwise we would have e = si = S2 which, by Lemma [2.111 would imply 
e = Si = S2 = — 1, obtaining some member of the family by Corollary 12. 131 Hence, 
we can suppose that e is different than the value described by formula (1261) above, and we 
conclude that T% 7^ 0. 

STEP #3: Clearly, one cannot expect to recover a unique / from e in general, as formula 
(125]) might suggests. Indeed, there are complex Hadamard matrices in which s\ = e, but 
S3 7^ /. The reason for this phenomenon is that formulas f l22|) and (j24D might be linearly 
dependent. After plugging fl23|) into ff24|) we obtain the expression 

FsGi - + (^3^2 - ^Gz) f = 0, 

which can lead us to one of the following three cases: 

CASE 1: Both the polynomials J-^Gi — F1G3 and J-3^2 — FiGz vanish identically, indepen- 
dently of e, meaning that in this case we do not have another condition on /. Luckily this 
can never happen, as by calculating a Grobner basis (again, to speed up the computations 
we have added the conjugates of the equations as well) we find that the polynomial 

(6 - l)b 2 (b - c)c(bc -d)(l + c + d)(b- d 2 )£{c, d) 

is a member of the basis. Therefore we have either c + d = — 1 or one of the degenerate cases 
described in the Canonical Transformation. The case c + d — — 1 implies that the set {c, d} 
consists of nontrivial cubic roots only. By directly solving the corresponding equations, we 
find that the quadruples (1, 1, to, u 2 ), (u 2 , u, u, u 2 ) and (1, 1, u 2 , u), (uj, u 2 , u 2 , lo) can vanish 
both polynomials, however these cases were excluded by the Canonical Transformation. 

CASE 2: Both the polynomials J-3^1 — F1G3 and J-3^2 — -^2^3 vanish for some |e| = 1, 
meaning that in this case we do not have another condition on /. However, having these 
numbers e at our disposal we can recover the two possible values of / from (1231 . Once we 
have the candidate pairs (e, f\) and (e, ^2) we readily calculate the remaining pairs (51,53) 
and (S2, S4) through the Decomposition Formula. Store all suitable sextics (e, si, S2, /, S3, S4) 
in the solution set SO Lb- Proceed to step #6. 

CASE 3: One cannot set (or have not set in Case 2) a unimodular e to make these two 
polynomial vanish at the same time. Hence we can derive formula fT231) for F(e). Proceed to 
step #4. 

STEP #4: Next we need to ensure that / is of modulus one. To do this, we calculate the 
fundamental polynomial 

V a ,b,c,d( e ) = 1-^3^1 — J^lGsf — \J 7 sG2 — ^2^3| 2 • 

After some calculations, it will be apparent that V has the following remarkable structure: 
V = a 8 b 8 c 6 d 6 P + a 8 b 8 c 6 d 6 P(2 + 2a + 2b)e + a 8 b 8 c 6 d 6 Qe 2 + Re 3 + Qe A + P{2 + 2a + 2b)e 5 + Pe 6 
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where the coefficients P, Q and R depend on the quadruple (a, b, c, d) only. If V = then 
the construction fails. Otherwise we find all possible roots of V[e) of modulus one. 

STEP #5: If the number on the right hand side of (1261) is not of modulus one, then the 
role of the pairs (e, /), (si, s 3 ) and (s 2 , s 4 ) is symmetric, and from all of the roots of V of 
modulus one we select all possible triplets (e,si,s 2 ) satisfying e + Si + s 2 = —1—a — b, which 
is needed to ensure orthogonality of the first two rows. From fl25|) we compute the unique 
companion values / = F(e), s 3 = F(s\) and S4 = F(s2). 

Otherwise, should the number on the right hand side of ( |26|) is of modulus one, then for 
every root e of V we calculate its unique companion value / = F(e), and then we use the 
Decomposition formula to determine the pairs (si, s 3 ) and (s 2 , s 4 ). 

If no unimodular roots are found at this point then the matrix E(a,b,c,d) cannot be 
embedded into a complex Hadamard matrix of order 6, and we do not proceed any further. 

At the end of this step we store all obtained sextics (e, s\, s 2 , /, S3, S4) in the solution set 
SO Lb- Typically two sextics are found. 

STEP #6: In this step one constructs the first three columns along the lines of steps #2-#5 
described above and obtain the set SOLq in a similar way. 

STEP #7: For every candidate solution from SO Lb and SOLq we check if the first three 
rows and columns are orthogonal, disregard those cases in which the submatrix B is singular 

and finally use Lemma to obtain the lower right submatrix D. Note that by Lemma \2. 141 

(3) 

we disregard members of the family K$ only. We check if D is composed of unimodular 
entries. 

STEP #8: Finally, we OUTPUT all unimodular matrices found during the process. We remark 
here that by Corollary 12.101 if no unimodular matrices were found then the submatrix E 
cannot be embedded into any complex Hadamard matrices of order 6. If unimodular matrices 
are found, then typically we find two matrices, as the solution set SOLb and SOLq contains 
two suitable sextics each, however, experimental results show that for each sextic in SOL B 
there is a unique sextic in SOLc making D unimodular as required. 

We have finished the discussion of Construction 13.11 The results are summarized in the 
following 

Theorem 4.1. Start from a submatrix E as in ^ and suppose that there are only finitely 
many (invertible) candidate submatrices B and C such that the first three rows and columns 
of the matrix G q (see §Z§) are orthogonal. Then Construction ^. 1\ gives an exhaustive list of 
all complex Hadamard matrices of order 6, up to equivalence, containing E as a submatrix. 

The interested reader might want to see an example of generic Hadamard matrices which 
can be described by closed analytic formulae, that is for which the fundamental polynomials 
'Pa,b,c,d an d V a ^b,d are both solvable. Such a matrix can be obtained when we choose the input 
quadruple (a, a, c, a) where the real part of a is the unique real solution of 49ft [a] 3 — 23? [a] + 1 = 
and c = (—a 3 + a 2 + a + 1)/ (a 4 + a 3 + a 2 — a). It is easily seen that the matrices we obtain 
starting from the submatrix E(a, a, c, a) are inequivalent from S^ and do not belong to the 
family 

Remark 4.2. When V = then the main difficulty we are facing with is that we have 
infinitely many candidate submatrices B. In this case we have the trivial restriction fl3]) on 
e, while the companion value / coming from (|2^|) is unimodular unconditionally. Although 
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in principle we can find three orthogonal rows through the Decomposition formula for every 
suitable e, we do not know which one to favourize in order to obtain a unimodular submatrix 
D via formula (1201) . Also, it might happen that the polynomial V a , c ,b,d obtained during step 
#6 shall vanish as well bringing another free parameter into the game making things even 
more complicated. In contrast, if both V a ,b,c,d ^ and V a , c ,b,d ^ 0, then we have a finitely 
many choices for the submatrices B and C and we can use Corollary 12.101 to conclude the 
construction. 

Remark 4.3. The polynomial V formally can vanish when we have T^Qi — = TzGi — 
TiGs = 0, however this is excluded by the Canonical Transformation and explained in details 
in Case 1 of step #3. It might vanish for some other, non-trivial quadruples as well making 
the whole construction process fail. In theory, the common roots of the coefficients of V can 
be calculated by means of Grobner bases, but as these coefficients are rather complicated 
obtaining such a basis turned out to be a task beyond our capabilities. Nevertheless, we 
conjecture that the case V = can be excluded completely in a similar fashion as we 
disregarded various quadruples during the Canonical Transformation. This would mean 
that all complex Hadamard matrices of order 6, except from Sg and K^ 3 \ can be recovered 
from Construction 13. ll 

It is reasonable to think that every complex Hadamard matrix of order 6 has some 3x3 
submatrix E leading to nonvanishing fundamental polynomials. In particular, we do not 
expect any complex Hadamard matrices of order 6 (except maybe and K^) which 
cannot be recovered from Construction 13.11 Therefore we formulate the following 

Conjecture 4.4. The list of complex Hadamard matrices of order 6 is as follows: the isolated 
matrix Sf\ the three-parameter degenerate family and the four-parameter generic family 
as described above. 

It would be nice to understand the structure of more thoroughly and express the 

entries of these matrices by some well-chosen trigonometric functions in a similar fashion as 

(3) 

Kg is described, however, as we have encountered sextic polynomials already the appearance 
of such formulas is somewhat unexpected. Also, it is natural to ask which matrices satisfy the 
conditions of Corollary 12. 181 An algebraic characterization of these matrices might lead to a 
deeper understanding of the generic family G^ and hopefully to the desired full classification 
of complex Hadamard matrices of order 6. 
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